Magnetic catalysis in hot and dense quark matter and quantum fluctuations 
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We analyze chiral symmetry breaking in quark matter in an external magnetic field at zero 
and finite temperature and quark chemical potential. We first give a brief overview of analytic 
results within the mean-field approximation. There the critical temperature for chiral restoration 
is increased by the magnetic field effect. Then we investigate the effects of matter and quantum 
fluctuations on the Magnetic Catalysis. More specifically, we compute the critical coupling as a 
function of the magnetic field and the temperature for zero and finite quark chemical potential in 
the presence of quantum fluctuations. As soon as a non-zero temperature and/or density is turned 
on, long-range correlations are screened and the critical coupling is no longer vanishing. We extend 
our dynamical results beyond the leading-order bubble resummation which results in a non-local 
four-Fermi coupling. This includes in-medium meson effects on the more quantitative level. 

PACS numbers: 11.30.Rd, 11. 10. Hi, ll.10.Wx, 21.65.Qr, 12.38.-t 



I. INTRODUCTION 

The state of matter in a strong magnetic field B still 
poses interesting theoretical as well as experimental chal- 
lenges. In the Landau quantization approach to these 
situations the energy dispersion relation of charged par- 
ticles is quantized along the directions perpendicular to 
B. There are a countless number of interesting phe- 
nomena related to the Landau quantization such as the 
Quantum Hall Effect in condensed matter physics [l[ , the 
Chiral Magnetic Effect in the relativistic heavy-ion col- 
lision 0, y] and related phenomena |J|, the primordial 
magnetic field in cosmological phase transitions [5]. In 
the present work we concentrate on the Magnetic Catal- 
ysis, [6], in the context of hot and dense quark matter. 

In relativistic heavy-ion collisions the state of mat- 
ter undergoes a phase transition or cross-over from the 
quark-gluon plasma of strongly interacting quark and 
gluon (quasi)-particles to the hadronic phase with chi- 
ral symmetry breaking and confinement. Recently, B- 
induced phenomena within these phase transitions are 
attracting more and more theoretical and experimental 
interest. This is motivated by the recognition of the pres- 
ence of extraordinary strong B in the early stage of non- 
central collisions [2] . The specific importance of B over 
the electric field E is that B is not screened by charged 
particles unlike E. The strength of B can get as large 
as \eB\ ~ Aq CD , while its life time is also characterized 

by ~ Aq^ d and thus is very small Q- Nevertheless, it 
is conceivable to anticipate observable B effects in the 
distribution of produced particles [8| since the produc- 
tion process occurs within the time scale of the strong 
interaction or even faster. 

By now model analyses according to the hydrody- 
namic description of the heavy-ion collision suggest 
that the thermalization is rapidly achieved before r ~ 
0.6 fm/c Q. Because B decays not exponentially but 



by the power as a function of time, the magnetic field 
strength is still comparable to the QCD scale when the 
system gets thermalized. Consequently it is a very inter- 
esting question how the QCD phase transition is modified 
by the strong B effects. 

Within the framework of effective model studies it was 
found that the chiral phase transition is significantly de- 
layed (i.e. the critical temperature is increased) by B, 
while the deconfinement transition or the Polyakov loop 
behavior is hardly affected |10l - [l2l |. This latter result on 
deconfinement, however, turned out to be inconsistent 
with the lattice-QCD simulation [lj, [l4|. Even though 
there are some attempts to account for the significant 
modification in the Polyakov loop behavior by introduc- 
ing extra couplings |15lll6l|. the microscopic details about 
the coupling between the gluonic dynamics and the exter- 
nal magnetic field are not fully understood yet. The miss- 
ing coupling is possibly attributed to screening by polar- 
ization effects, which has the same origin as the failure of 
chiral model approaches to baryon-rich matter [13, [l8| . 
There it has been shown within QCD computations that 
the polarization effects have a considerable impact on the 
location of the confinement-deconfinement phase transi- 
tion; see Refs. [19|, [20] ■ Indeed, in Ref. [2l[ it has been 
shown that the interplay between strong electromagnetic 
and chromo-electric and chromo-magnetic fields gives rise 
to polarization effects that catalyze the deconfinement 
phase transition. This already hints at the similarity be- 
tween large- £? effect and high-/i q effect which may give 
a clue for attacking the problems at high density, as em- 
phasized in Ref. |22| . 

In contrast to deconfinement, the B effects on the chi- 
ral phase transition is naturally understood from the 
qualitative effects implied by the so-called Magnetic 
Catalysis 6] which was first discovered in fermionic the- 
ory with four-Fermi interaction, and later the argument 
was extended to QED as well [23] . When B is sufficiently 
strong, the magnetic field plays a role as the catalysis 



that induces a finite chiral condensate (ipi/>). Since chi- 
ral symmetry breaking is enhanced by B, it is natural to 
anticipate a higher critical temperature for chiral restora- 
tion than in the case with vanishing B. This expecta- 
tion was indeed confirmed in recent model studies [24] 
and indeed was observed already before the discovery of 
the Ma gnet ic Catalysis |25|. In recent lattice-QCD sim- 
ulation [14[ it was found that T c significantly decreases 
for larger B, which can be interpreted as a result of the 
entanglement with the Polyakov loop that is affected by 
the magnetic screening [26j. Such screening has been 
discussed in Ref. [21( for the gluonic potential. In this 
sense, as long as only the chiral sector is concerned, the 
enhancement of the chiral condensate by B is not contra- 
dictory to Ref. [14]. Also, regarding another possibility 
of the inverse Magnetic Catalysis at finite density, see 
Ref. H3- 

The purpose of this paper is to study the B effects on 
the chiral phase transition including quantum fluctua- 
tions. These effects of quantum fluctuations are quite 
naturally included with renormalization group (RG) 
flows which also shed light on the physics mechanisms 
at work (see Refs. [l2|, [2a [29j for related works). We 
calculate the /3-function for the four-Fermi coupling con- 
stant Afe as a function of an IR cutoff scale k. At large 
cutoff scales k implying large energies we are in the chi- 
rally symmetric phase with vanishing chiral condensate 
(ipip) = 0. Chiral symmetry breaking is signaled by a 
singularity in the four-Fermi coupling related to the ap- 
pearance of massless pion (and sigma) propagation in the 
intermediate state of the four-point vertex at vanishing 
momenta. Hence by decreasing the cutoff scale k the 
RG-flow of Afe hits this singularity at some point (see 
Ref. [30] for a recent comprehensive review). In this pa- 
per we adopt a chiral model for simplicity. In this way we 
can get a simple formula for the critical surface in space 
spanned by the four-Fermi coupling Xk, the temperature 
T, and the quark chemical potential fi q with and without 
a constant magnetic field B. The extension or embed- 
ding of the model results obtained here within the RG 
approach to QCD is straightforward; see Refs. [20ll30j. 

This paper is organized as follows. We first present 
a brief overview of the mean-field results for the criti- 
cal coupling in Sec. [TTJ This part is used for introducing 
our notation and it also has the benefit of allowing an 
easy comparison of the mean-field arguments and results 
with those of the RG analysis in Sec. IIIII Finally we ex- 
tend the RG analysis by including meson effects beyond 
the leading-order bubble resummation to the four-Fermi 
coupling in Sec. IIVI As we shall see in Sec. IIIII there 
are three different kinds of Feynman diagrams that con- 
tribute to the /3-function of Xk- Two of them tend to 
break chiral symmetry and are genuinely fermionic in na- 
ture. The last one forms a box-type diagram that favors 
chiral restoration and this box diagram contains interme- 
diate meson states with finite momentum insertion. This 
additional resummation establishes the well-known effect 
that fermionic fluctuations tend to break chiral symme- 



try whereas mesonic ones tend to restore it. We find that 
the box diagram has only a minor effect on the Magnetic 
Catalysis when B is large. Section |V| is devoted to the 
conclusion and outlook. 



II. OVERVIEW OF THE MEAN-FIELD 
CALCULATION 



Here we recall spontaneous chiral symmetry breaking 
on the mean- field level with and without magnetic field. 
We intend to capture the generic feature of chiral sym- 
metry breaking in an effective description. Integrating- 
out the gluons in QCD as well as the quark fluctuations 
with momenta larger than a ultraviolet (UV) cutoff scale 
A would lead to a chiral effective theory formulated in 
terms of quarks with all kinds of interaction vertices; see 
e.g. Refs. pjl 120, biH |32[ . Among these vertices only the 
point-like four-point interaction is retained for simplic- 
ity, which defines the Nambu-Jona-Lasinio (NJL) model. 
In what follows we consider hot and dense quark mat- 
ter with three colors, N c — 3, and two massless flavors, 
JVjf = 2. Then the Lagrangian density at the UV scale A 
reads 



c K = $W> + \$l 



"Vf r A $£ #*V2f 



(i) 



where a and j3 run in color space, a, b, c, d in flavor 
space, and i, j, I, m refer to the Dirac indices. We take 
the sum over repeated indices as usual. For the concrete 
form of ^Aiffm we assume an equal mixing of the U(1)a- 
symmetric and U(l)A-breaking terms so that the iso- 
scalar pseudo-scalar particle (770) and the iso- vector scalar 
particles (do) decouple from low- lying spectra. The re- 
maining interaction vertices contain the iso-scalar scalar 
channel (ex) and the iso- vector pseudo-scalar (7?) channels 
only, 

Ta^ = X A [S ij S lm S ab S cd + (h 5 h(h 5 )i m (r n ) ab (r n ) cd ] , 

where r"'s are the Pauli matrices in flavor space and n 
runs from 1 to Nf — 1 corresponding to pion degrees of 
freedom. 

We close with the remark that the Lagrangian ^ was 
introduced as already comprising all quantum effects of 
quark fluctuations with momenta larger than the cutoff 
scale as well as all gluonic fluctuations. This implies that 
(J2|) should have included terms depending on A as well 
as the physical mass scales of QCD (Aqcd, ^quark) such 
that if introducing the missing quantum fluctuations of 
the quarks, the result is A-independent. This necessity is 
avoided by assuming A to be such that these terms are 
minimized. This gives a physical meaning to A and it is 
rather Aqcd up to a constant than a UV cutoff scale. 



A. Gap equation and the critical coupling constant 

In the mean-field approximation the four-Fermi inter- 
action is decomposed into the mean-field (tptp) and the 
fluctuations about it. We pick up only the contribution 
from the non-crossing color structure, (■ip a ip a ) -tp^ip 13 , but 
simply neglect the crossing one like (■ip a ipP) -ip^ip a . Such 
a treatment is justified in the large- N c limit in which the 
mean-field approximation also becomes exact. The con- 
stituent quark mass, M — —\\(ijjip), appears as a result 
of the spontaneous chiral symmetry breaking, and the 
mean-field Lagrangian density reads 

Aa 



£mf = iptfhp + Aa(V'V') ipip — 7r{tyi>) 



ip(i@-M)ip 



M 2 
2Al 



(3) 



in the chiral limit. The Lagrangian density ([3]) bilinear 
in the quark fields, and hence the path integral quanti- 
zation is done by a Gaussian integration. The resulting 
thermodynamic potential is a function of M and is given 
by 



ptt/V = -2N C N { f t^Uw + Tln[l- 

J p 2< A 2 (2tt) j I 

+ T\n[l + e- l3{uj+ ^]\ + 



(w-/j q )n 



M 2 
2Al 



with the quasi-particle energy dispersion relation, uj = 
^Jp 2 + M 2 . The momentum integral is cut at the cutoff 
scale as the mean field Lagrangian © supposedly aris- 
ing from integrating-out the quarks with momenta larger 
than the cutoff scale, p 2 > A 2 . The overall constant 2 
comes from the spin degeneracy. For the moment we 
restrict ourselves to vanishing temperature and density; 
T = Hq = 0. Later we extend the analysis to include 
thermal and medium effects. In this simple setup the 
momentum integral in Eq. (U) can only lead to a func- 
tion of £ = M/A only, up to a dimensional factor A 4 . The 
zero-point energy contribution gives the leading contri- 
bution to the thermodynamic potential in Eq. Q and 
can be evaluated explicitly as, e.g. |33| . 

d 3 p 



2N C N { 

)p 2 <A 

N c N f A 4 
8tt 2 

NcNfA 4 



2tt 



v / T+^(2 + £ 2 



e 



In 



v^+e 



Vi+F+i 



4tt 2 



i+e)+o(e) 



(5) 



Note that we do not loose the generality of our analysis 
by assuming small £ = M/A because we are only inter- 
ested in the onset of chiral phase transition where £ starts 
taking a finite value from zero. With the expansion in £ 
the thermodynamic potential simplifies, to wit, 



pn/v 



N C N { A 4 
4tt 2 



1 



1 - 



2k 2 



N C N { A A A 2 






(6) 



where we have dropped the higher order terms in £. The 
critical coupling associated with the second-order phase 
transition should correspond to a point at which the cur- 
vature of the thermodynamic potential crosses zero. In 
order to have a non-zero £ the potential curvature should 
be negative and this condition leads to an inequality [33j , 



A A A 2 > 



2tt 2 



(7) 



This is the condition for the NJL model to accommodate 
for spontaneous chiral symmetry breaking. Seemingly 
this depends on the UV cutoff scale A. However, as we 
have already argued above, A relates to Aqcd in the 
current effective field theory set-up. 

Finite temperature effects are taken into account in a 
straightforward extension. The full dependence includ- 
ing the finite-/i q effect will be considered in the numerical 
analysis in Sec. IIIII Again we aim at simplicity and ana- 
lytically estimate the finite-T corrections. This leads to 
simple understanding of how the chiral phase transition 
takes place within this effective model. Since we are in- 
terested in the regime with infinitesimal £, we can make 
use of the high-T expansion for T ^> M, which yields the 
following thermodynamic potential, 



(4) pn/v 



N C N { A 4 
4T 2 



2^ 2 



N C N { A A A 2 



TV" 

y 



(8) 
Then, in the same way as the condition (|7|) has been de- 
rived, we see that spontaneous chiral symmetry breaking 
requires, 



A A A 2 > 



2tt 2 

AUVf 



tI?i ?2 



(9) 



In the effective model studies, usually, the four-Fermi 
coupling is fixed by observables in the vacuum. For a 
given Aa that satisfies the above inequality, we can read 
the critical temperature directly from Eq. ([§]) as 



T r 



'3A 2 



6 



N c N f \ A 



(10) 



We readily confirm that Eq. ([10]) gives a conventional 
value for the critical temperature, i.e. T c = 175 MeV 
by plugging-in the standard NJL-model parameters, 
e.g. [lif: A = 631 MeV and A A /2 = 0.214 fm 2 . 



B. Effects of the magnetic field 

We proceed by introducing an external magnetic field 
B. In the present work we restrict ourselves to spatially 
homogeneous and temporally constant B fields. In the 
presence of B (along the z-axis) the transverse momenta 
(i.e. the x and y components) of spin-1/2 particles with 
electric charge q are quantized into the Landau levels and 



the momentum integration with spin sum of a general 
function f(u>) is replaced as 



d 3 p 



/(«) 



(2tt) 3 

_\ q B\;^ f di>. 



j)2<A2 v^'T JA,S 

N A 



d 3 p 



/(w„ 



(11) 



where 2 in the left-hand side is the spin degeneracy, which 
is replaced by the spin-degeneracy factor a n that is 1 
for n — (Landau zero-mode) and 2 otherwise. In the 
second line, p 2 = p 2 z +2\qB\n and N AiB = 9 g [A 2 /(2\qB\)], 
where 9 g [n + x] = n for n £ Z and < x < 1 . 

Now, let us continue our analytical consideration tak- 
ing the simplest limit of large magnetic field, more specif- 
ically 2\qB\ > A 2 . In this particular limit only the Lan- 
dau zero-mode contributes as 8 g [A 2 /(2\qB\)} = and 
the lowest Landau level approximation (LLLA) is exact. 
Note that this should be already a rather quantitative ap- 
proximation for heavy-ion collisions as the magnetic field 
strength there is of order of Aqcd which is related to our 
cutoff scale A. In any case we shall also consider general 
magnetic field within our numerical calculations. In the 
LLLA, at /i q = 0, we can express the thermodynamic 
potential as 



w-*^/! 



f=u,d 



Wo 



2Tln(l + c 



-/3u> \ 



M 2 

2Aa" 



(12) 



Here we defined luq = \/p 2 + M 2 and the quark electric 
charges are q u = (2/3)e and qd = — (l/3)e. From this 
form of the thermodynamic potential one can easily un- 
derstand the essence of the Magnetic Catalysis. At T = 
we can expand the zero-point energy for small £ = M/A, 
as we did previously, to find the thermodynamic poten- 
tial up to 0(£ 2 ) as 



pn/v 



-JV, 



E 

f=u,d 



\lfB\A 2 
4tt 2 



^+-Vr +^e 



2A, 



(13) 

As we have already argued, the sign of the curvature de- 
termines whether chiral symmetry is spontaneously bro- 
ken or not. The remarkable feature in the LLLA is that 
ln(l/£) in the coefficient of £ 2 is positive and arbitrar- 
ily large for £ — > and thus it can eventually overcome 
A 2 /(2Aa) from the last term. This means that chiral 
symmetry is always broken however small A a is. The im- 
portant point is that such a logarithmic term appeared in 
the coefficient of £ 4 in Eq. ([5]) when B was not applied. 
In the LLLA with strong B the phase-space volume is 
suppressed and £ 4 is reduced to (|q/-B|/A 2 )£ 2 . Because 
ln(l/£) changes its sign from positive to negative for large 
£, the thermodynamic potential (I13|) has a minimum as 



a function of £ that solves the gap equation. In fact, by 
taking the derivative, we can write down the gap equa- 
tion, -A^ c X; / (|g/S|A 2 /27r 2 )ln(2/£) + A 2 /AA = 0, whose 
solution £o or the corresponding constituent mass is 



M = £ A = 2Aexp 



2tt 2 



Nc\AEf\VfB\ 



(14) 



This expression for the condensate is quite analogous to 
the superconducting gap in the BCS theory. We can see 
that the critical temperature at which the condensate 
melts is also given by the same relation as the one known 
in the BCS theory. 

We can carry out the high-T expansion in the same 
way as previously using 



J-oo 2tt V 



-fiu> \ 



irT A 2 /, A£ 

1 In — - + 7 - 

12 47TTV ttT ' 



up to 0(£ 2 ). Interestingly enough, a logarithmic term 
proportional to £ 2 ln(l/£) appears from this finite-T in- 
tegration as well as the vacuum contribution that we al- 
ready mentioned. We can make sure that these loga- 
rithmic singularities exactly cancel out and the Magnetic 
Catalysis is lost at finite T. This trend had been reported 
in the literature [25| and was found also in the context 
of the Magnetic Catalysis [3J]. In the RG analysis in 
Sec. IIIII the cancellation of the IR singularities will be 
reconfirmed in a more transparent manner. 

The curvature of the thermodynamic potential no 
longer has a logarithmic singularity at finite T and the 
critical coupling constant becomes finite again. This is 
actually the reason why we can expect a chiral phase 
transition at a certain temperature even in the LLLA 
calculation with strong B. After all, the thermodynamic 
potential (fT5j) receives the finitc-T corrections as 



pn/v ■ 



/ 



\q f B\A 2 



In 



4tt 2 
/2e^A 



1 



vr 2 /T x 2 



3 \A 

2ir 2 



(16) 



V nT J N c X A J2 f \q f B 



e 



The zero of the coefficient of £ 2 gives a condition for the 
critical coupling or the critical temperature as obtained 
in Eq. (flOl) . In this case of the LLLA we arrive at 



T 2e^A 
l c = exp 



2^ 2 



^cA A E f l<Z/£| 



= —M 

IT 



(17) 



This relation between the condensate at T = and the 
melting temperature is perfectly identical to the results 
from the BCS theory. If the coupling constant Aa is 
infinitely large, all expressions appear very simple; the 
critical temperature without B is read from Eq. (p~0|) as 
T c = (V3/n)A ~ 0.551A, and that at strong B is from 
Eq. 02} as T c = (2e 7 /7r)A ~ 1.13A. Therefore the crit- 
ical temperature is indeed raised by the B effect, which 
is in qualitative agreement with effective model calcula- 
tions. 



III. RENORMALIZATION GROUP ANALYSIS 

In the previous section we have discussed the stan- 
dard mean-field analysis where quark quantum fluctua- 
tions are taken into account at one loop. In the present 
section we extend the analysis beyond the one-loop level 
by means of RG-techniques. This also allows us to for- 
mulate the physics mechanism of the chiral phase tran- 
sition from a different and rather simple point of view. 
To that end we here introduce an IR-cutoff term into 
the Lagrangian and derive a functional RG equation for 
the free energy or effective action, the Wetterich equa- 
tion 35]. In this approach the kinetic term is modified in 
momentum space in order to suppress momentum modes 
below the given infrared (IR) cutoff scale k. This leads 
to an infrared regularized free energy or effective action 
Tk[ip,ip] which reduces to the thermodynamic potential 
fi at vanishing IR cutoff k = on the equations of mo- 
tion (EoM) of ipj'ijj, that is fl — T/^oIeoM- For the sake 
of simplicity we resort to a three-dimensional IR-cutoff 
which allows us to perform the Matsubara summation 
analytically. Then the fermion propagator in Euclidean 
space takes the following form, 



G(p) 



1 



74P4 + 7-PI 1 + rk(p)} 



(18) 



Here r k {p) is the cutoff function with r k (p) ~ k/\p\ for 
p 2 < k 2 and r k (p) — > for p 2 > k 2 . Note also that 
the dimensionless shape of r k is genuinely a function of 
|p|/fc. Since we are interested in the chiral properties of 
QCD, we have avoided implementing the cutoff function 
as a momentum-dependent mass term but in the kinetic 
term. The former cutoff explicitly breaks chiral symme- 
try and would have similar effects on chiral properties 
as the Wilson mass term has on the lattice. The scale- 
dependence of Tk is then encoded in the flow equation, 
schematically written as 



7 -pdtTk 



a t r fc [^,v] = -Tr[-^ T) 



I\. ' '[ip,ip] +j-pr k 



^(1,1) 



V 



n ' WH^r^fe, 



64> 



'5fl> 



(19) 



where t — ln(fc/A) is the logarithmic IR scale with ref- 
erence scale A and the trace Tr sums over momenta, 
Dirac indices, flavors. (For QCD-related reviews, see 
Refs. [30, [3a - l39j V In the present work we ignore the 
effects of the multi-scattering of quarks, which is phase- 
space suppressed. We rush to add that this argument 
gets weaker at higher densities. For the time being we 
make the additional approximation of a local effective 
four- Fermi vertex X k , which will be relaxed in the next 
section. This leads us to the following simple Ansatz for 
the effective action, 



r fe = di i/n< 



+ ^r^ a nt^f^ d J) (->()! 



where ^kijfm is given by Eq. ([2]) with A — !• fc. This form 
is substituted into the flow equation (fH?|) and we take 
the functional derivatives. The left-hand side of the flow 
equation is then, 



d t if < 2 > 



d t 



7 7 



-IV 



V 



V 



= 5{x — y) S(y — z) S(z — w) 

x (S^S^d t T k t^i - S aS S^d t T k ^) . (21) 

The right-hand side is a bit complicated. There arise 
twelve terms, six of which are proportional to 8 a ^8 1 and 
other six proportional to S^S^ 1 , respectively. Here, let 
us focus on the part of 8 a ^8 lS only. In addition, because 
we aim to have the flow equation for A&, we can sim- 
plify the calculation by contracting the Dirac and the fla- 
vor indices of external legs, namely, SijSi m 5 ab 5 cd T kfj^ = 
16N 2 X k . Then six terms are no longer independent but 
only three remain distinguishable. The flow equation is 
eventually given by 

d t T k ?Sr = 2N c Tv[G ]m T b ^ l G pq Tf r ^G rs (i-p}s 1 d t r k ] 

— 2Tv[Gj Tn r 7 ^° ip G pq r qr '[[G rs (^ ■ p) S j d t r k \ 

- 2T r [G 3m r b ^G pq r d q ff r b G rs (l ■ P)s 3 d t r k ] , (22) 

when no momentum is carried by any external legs. In 
Eq. (f2"2")l we have also dropped the subscript k on the 
vertices. There are four different types of Feynman dia- 
grams that contribute to the flow equation as displayed 
in Fig. Q] Here (a) corresponds to the first term in the 
right-hand side of Eq. (|2"2"]) . (b) and (c) to the second 
term, and (d) to the third term. 

After some calculations we can arrive at the following 
form of the flow equation with N{ = 2; 



d t \k 




\p\d t r k (p) 



(23) 



where the factor four in Eq. (|2U|) originates in the trace of 
the Dirac index and the factor 2N C comes from (a) which 
is the leading diagram in the 1/N C expansion. The factor 
3/4 = 1 — 1/4 is a combination of 1 from the diagrams 
(b) and (c) and 1/4 from the diagram (d). 

Before we come to the explicit computations we would 
like to discuss the structure of the flow in Eq. (|2"3")l as 
it already reveals a lot of interesting properties of the 
flow and is also important for the computations and the 
interpretation of the results. The flow (J23I) has a total 
^-derivative structure and hence leads to results at k = 
that are independent of the chosen regulator r k up to (re- 
)normalizations at the initial scale A. To see this more 
clearly we rewrite it as a flow for 1/Afe, to wit, 



di 



A/,. 



-9 t [n fc (o)-n A ( 



(24) 




FIG. 1. Diagrammatic representation of the flow equation for the four-Fermi interaction. At each vertex the color index is 
contracted along the connected lines. The ellipses are neglected contributions from higher loop diagrams beyond the Ansatz (|20|) 



with 

n fc (o) 




(25) 



In Eq. (|24[) we have introduced the subtraction term for 
having explicitly finite expressions before ^-derivatives 
are taken. In its form (|24[) the flow is easily integrated 
and leads to 



i- = -^-[n fc (o)-n A (o)] 



(26) 



We immediately deduce that the dependence on the reg- 
ulator drops out at k = for consistent choice of Aa; 
evidently there is no dependence on r (p) = 0. In turn, 
1/Aa + IIa(0) is required to be A- independent, otherwise 
Ao would depend on A. The latter condition encodes 
RG-invariance of the theory, not to be mixed-up with 
renormalizability. We also deduce from Eq. (|26"|) that in 
the leading order 1/Aa oc — IIa(O) = c(rA)A 2 . The pref- 
actor depends on the chosen regulator and hence defines 
the regularization scheme. This is seen more clearly if 
comparing the dimensionless coupling, 



At — Ai-A; 



at the initial scale A for two different regulator 



(27) 

„M J 2 ) 



This leads to X\ 



>/Af 



= A 2 ) /M) 



/c' K For agreeing physics 



scales the above ratio should be one and we identify 
A< 2 ) = c (1) /c^A with A (1) = A. In other words, we ex- 
pect that the equivalence of the mean-field inequality 
depends on the chosen regulator simply by a multiplica- 
tive factor in front of A; 



A A A 2 > 



c(r A ) 2tt 2 

c(r shar P) NcNf ■ 



(28) 



where c(r s A a ' lp ) stands for the implicit regulator used 
in the mean-field computation, that is the three- 
dimensional sharp cutoff; rl arp (p 2 ) = 9(p 2 — fc 2 ). We 
emphasize once more that this only reflects the depen- 
dence of Aa on the chosen RG-scheme. 

In summary we conclude that Ao does not depend 
on r\ either and any regulator-dependence is removed. 
We emphasize that this is a particularity of the present 
approximation and hinges on its explicit total deriva- 
tive structure. Such approximations are singled-out as 



regulator-independent and also related to optimization 
criteria that amongst other properties (re-)enforce total 
derivative structures for the flows; for more detailed dis- 
cussions, see Ref. [38, 40]. For later purpose we also 
resolve Eq. (|2l>| for Afc, 



Aa 



i-A A [n fc (o)-n A (o)] 



(29) 



Equation (f2l)|) elucidates the resummation structure of 
the flow. It encodes a bubble resummation of the di- 
agrams in Fig. [TJ This is seen within a diagrammatic 
expansion of Eq. (|2"6")l in orders of Aa, 



oc + )oo 



(30) 



If we interpret the dashed line on the left-hand side as 
a meson propagator it results simply from the bubble 
resummation shown in Fig. 1301 The first term on the 
right-hand side, i.e. the tree-level diagram, gives Aa and 
hence a tree- level meson propagator Aa ■ The second term 
comes from the one-loop polarization, that is, Xjjlk(0), 
where the momentum argument is zero because we only 
consider the situation that all external legs bring no mo- 
mentum insertion. If we lower the temperature towards 
the hadronic phase the four- Fermi vertex function has to 
diverge at the phase transition as it contains a massless 
sigma and pion propagation. Indeed, when introducing 
effective mesonic degrees of freedom with the a and -k 
coupling to the related quark bilinears, it can be shown 
that the meson mass is proportional to 1/Afe. This makes 
it clear that a single pole in the four-Fermi coupling does 
not merely signal the validity bound of the approxima- 
tion but rather the onset of condensation. 

This ends our discussion of the formal properties of 
the flow ((23)) and we proceed on how to practically solve 
it. As we have argued, its solution does not depend on 
the chosen regulator. This fully justifies the choice taken 
already in restricting ourselves to three-dimensional regu- 
lators. It allows us to perform the p^ integration or sum- 
mation with respect to the Matsubara frequency most 
easily. The result is 



TV - = ^- + ±J— 

Z-* {p\ + u 2 ) 2 Auj 2 duj 2w e^ 



1 



(31) 



with u = \p\\l + rk{p)\- We still have to specify the 
three-dimensional regulator, rj,(p), to proceed further. 
We again utilize the regulator- independence and adopt 
the three-dimensional flat cutoff, [41], for the sake of sim- 
plicity, 



^(p) 



\P\ 



1 0(k/\p\ - 1) 



(32) 



Such flat cutoffs were first introduced in Ref. [42j for the 
purpose of optimization are also most useful for analyt- 
ical studies. It is called "flat" because it renders the 
spatial dispersion momentum-independent for momenta 
below the cutoff scale, 



(33) 



w(rjrw*/id - 1) = k 



where u> denotes a quantity defined just blow Eq. (1311) . 
For momenta \p\ < k, hence, the Matsubara sum in 
Eq. JSIJ turns into 



1 

2fc2 



1 



1 



/3k- 



■Jik 



3 /3fc + l " ( e /3fc + l)2 



(34) 



which is completely momentum-independent. The inte- 
grand in the four-Fermi flow (j2"3")l is also proportional 
to the scale derivative of the regulator, <9 t rjJ at (p) = 
(k/\p\)9(k/\p\ — 1) which limits the integration to mo- 
menta |p| < k where the integration is trivial. This leads 
us, particularly when B = 0, to 



ftn fc (o) 



2N r 



)N t 



3tt 2 



1 



c /3k + I 



/3k 



.,dA- 



(e 



pk 



I) 2 . 
(35) 
to be inserted in Eq. (|24p . In later discussions it also 
turns out that the choice (133[) for rk(p) avoids unphysi- 
cal oscillatory behavior in the flow from the discrete Lan- 
dau levels. For going beyond the present approximation 
which we shall do in the last section it is worth notic- 
ing that this cutoff is optimized in the high tempera- 
ture limit, [33, H2], and keeps some of this property also 
at momentum-independent approximations. For fully 
momentum-dependent approximations as well as finite 
magnetic field it introduces momentum non-localities 
(momentum flows in diagrams) which affect the quan- 
titative accuracy (see Ref. j43|). In short, within the 
approximations used in the present work the regulator 
works sufficiently well and has the advantage of analyti- 
cal tractability. Within extensions of this work we shall 
aim at quantitative precision and will resort to regula- 
tors which maintain momentum locality and minimize 
the momentum flow. 



A. Flow and chiral symmetry breaking 

We first concentrate on the relation between the flow 
of the four- Fermi coupling and chiral symmetry breaking 
at T — and B = 0. This gives access to the simple 



structure and understanding of chiral symmetry break- 
ing within the flow equation approach, and hence serves 
as a useful example for understanding the results in the 
presence of medium effects. In order to facilitate the di- 
rect comparison with the mean-field results, we shall first 
discuss the large-iV c limit at leading order. Then, only 
the diagram (a) in Fig. Q] contributes. Note that due to 
the simple full ./V c -dependence of the flow (see Eq. ([23|) ). 
the limit is easily undone by replacing N c — > N c + | in 
the large N c result. For the flat regulator and T = the 
flow of the vacuum polarization, Eq. (|3"5|) . is iV c /(37r 2 ) 
and we arrive at 



2 1.2 



<9 t A fe = - 



as well as 



N c Nf\ jk 
3tt 2 



Aa 



NcN f X A fh2 vJi 

6 9 



-(fc 2 -A 2 ) 



(36) 



(37) 



with a positive initial coupling Aa- From Eqs. (|36[) 
and (|37p we infer that lowering k increases A^. At 
N C N{ AaA 2 /(67t 2 ) = 1 the coupling diverges and has a 
single pole. We have already argued that Eq. (f2"5| has 
to be understood as the onset of spontaneous chiral sym- 
metry breaking, and the occurance of massless mesonic 
modes. Hence, with the present regulator the condition 
for spontaneous chiral symmetry breaking, Eq. 
given by 



A A A 2 > 



6tt 2 

N c N f 



(38) 



This result is the counterpart of the condition (|7|) in the 
mean-field theory. The discrepancy by a factor 3 relates 
to the different regularization scheme between the sharp 
cutoff in the mean-field treatment and the optimized cut- 
off in the RG analysis. We can identify, roughly speaking, 
the RG-scale A 2 in the present RG study with 3A 2 in the 
mean-field approximation. 

The same conclusion can be drawn from a fixed-point 
argument, directly utilizing Eq. (|36j) : by changing to the 
dimensionless coupling A in Eq. (|27| the flow (l36l) changes 
to 



p- x = d t \k = 2A fe - -jj^a" A fc 



(39) 



The /3-function in Eq. ([29) vanishes at A* = 67r 2 /(7V c A^ f ) 
which is an attractive UV fixed point. Of course, this 
defines exactly the critical coupling constant given in 
Eq. ([3"8")l . If Aa starts from a value above A*, the flow 
grows larger for smaller k, which indicates the sponta- 
neous breaking of chiral symmetry 

With this simple picture of chiral symmetry breaking 
we proceed to finite temperature T and quark chemical 
potential /i q . More specifically we search for the critical 
temperature (or chemical potential) above which chiral 
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FIG. 2. Flow pattern of the four-Fermi coupling constant 
without B (upper figure) and with B (lower figure). In the 
latter case chiral symmetry is always broken as a result of 
dimensional reduction. 



symmetry is restored for a given Aa- The flow equation 
now receives medium corrections, to wit, 



<9*Afc — 2Afe — 



N c Nj\l 
3tt 2 



x \ 1 - (1 - d t ) [ra(fc-/i q ) + n(k+fi q )] 



(40) 



where n(k) = l/(e^ k + 1) is the Fermi-Dirac distribution 
function. Assume now that the initial coupling exceeds 
the fixed point coupling at vanishing temperature, i.e., 
Aa > A* with A* in Eq. (pJ5)) and hence (i^\k=K < 0. 
we again make contact to the mean-field analysis within 
a solution of the differential equation (l4Tfl) at /i q = 0. 
We can readily integrate it to write down the running 
coupling constant, 



A 



Aa 



1 + 



A^VfAA 
3tt 2 



V e^+l (e^'+l) 2 

(41) 
Here, because we are interested in T c inferred from the 
Landau pole associated with the flow of k — > 0, we make 
an approximation of k ~ and A ~ oo in the matter part 
and evaluate the integration as 



6k' k 1 



2/3k'e 



/3k' 



3 0k> + i ( e /3fe< + iy 



n 2 T 2 



(42) 



Then, the running coupling (|4i"j) can be well approxi- 
mated as 



A fc 



A A 



1 + 



N C N { A A (h2 x2 ^ r ,, 
6tt 2 



■(fc 2 -A 2 +^ 2 T 2 



(43) 



from which the Landau pole position leads to the condi- 
tion for the critical temperature; 



T r 



(i 



/A 2 



ir 2 N c N f A A 



(44) 



With the cutoff replaced as A 2 — > 3A 2 , this is precisely 
consistent with the mean-field result in Eq. (QI 



B. Magnetic Catalysis 

We introduce finite B using Eq. (111]) and inserting ap- 
propriate spin projection operators in the trace of the 
Dirac indices. In the LLLA we can easily understand the 
Magnetic Catalysis from the flow pattern. At T = the 
flow equation in the strong B limit takes the following 
form: 



N, 



^ = -;dr5>/^ 



(45) 



/ 



with t = ln(fc/A), which we can solve immediately as 

A A 



At- = 



i+^Ew 



(46) 



Because t = ln(fc/A) ranges from (fc = A) to — oo 
(k = 0), the denominator of Eq. (|46p inevitably hits the 
Landau pole in the flow of k in sharp contrast to the 
zero-.B case. The existence of the Landau pole implies 
that the chiral symmetric state is always unstable and so 
chiral symmetry should be broken in the LLLA. This is 
nothing but the manifestation of the Magnetic Catalysis. 
In the same way as previously, we can understand this 
from the behavior of the /3-function of A&. In contrast 
to Eq. (J39]) the mass dimension of Afe is now balanced 
by not k 2 but \q_fB\ because of the dimensional reduc- 
tion. Therefore, the first term in Eq. (|39l) . 2Afc, does not 
appear, and the behavior of the /3-function changes as 
sketched in Fig. [21 More specifically, introducing a di- 
mensionless coupling hy \ k = (J2t |9/-S|/27r 2 )Afe, the 
flow equation reads, 



J \m 



-N C (X 



(B)> 

* i 



(47) 



It is obvious from this that the flow goes to infinity 
regardless of the initial point of Aa- From the so- 
lution (|4"6")l we can easily locate the Landau pole at 



ln(fc /A) 



-2tt 2 /(N c Xa Ylf \QfB\), that characterizes 



the typical scale for the condensate, i.e. 



Mq oc fcp = A cxp — 



2-k 2 



NcXaE f \lfB\ 



(48) 



which gives an estimate for the chiral condensate up to an 
overall coefficient. We note that the exponential factor is 
identical with that in Eq. ([i"4")l since it does not depend on 
A which carries the scheme-dependence. It does depend 
on J2f \QfB\X\ which is already dimensionless. 



C. Finite temperature and density 

It is an interesting question how the chiral phase tran- 
sition at finite T and/or /i q is possible even in the LLLA 
in the RG formulation. In other words, the question is 
how the Magnetic Catalysis is lost and chiral restoration 
becomes possible when a finite T and/or /i q is turned on. 
The running coupling constant has extra contributions 
at finite T as 



A& = Aa 



1- 



JVcAa 
2tt 2 



x 1 



2 2f3k' e^' 

cf> k ' + 1 ~ (e^' + If 



(49) 



As we did before, we can approximate k ~ in the in- 
tegration range as long as we are interested in T c only, 
and also take A ~ oo, which makes the last term in the 
integration part in Eq. (|4"9")l as 



dfc' 



2(3k'e 



Pk' 



( e /9k' 



= 1 



(50) 



The first and the second terms are more interesting than 
the third one. The Magnetic Catalysis stems from the 
IR singularity from J, dk'/k' = — ln(fc/A), while the 
finitc-T contribution exactly cancels the IR singularity 
as explicitly checked as 



dfc' 



3)8*' 



1 



/3A da; c x - 1 
x e x + 1 



ln(1.13/3A). 

(51) 

This recovers the same expression for T c as Eq. (|17p up 
to the prefactor that again depends on the regularization 
scheme; from the zero of the denominator in Eq. (|49j) we 
can infer in the LLLA, 



T c = 0.42Aexp 



2tt 2 



AUa £/ \q f B 



(52) 



We can relax the LLLA and solve A^ for arbitrary B at 
finite T and /x q numerically. We can write the running 



coupling constant down as follows 
N C X A 



A/c = Aa 



1 



2tt 2 



Ei^i 



A dfc' 



E a » V 1 " ^^l 1 - n(k'-^) n{k'+^) 



+ kdk [n(fc'-/i q ) + n(k'+ii q )] j 



-i -l 



(53) 



In this expression n can take a positive integer up to 
Nk,B — @g[k 2 /{2\qfB\)]. From the zero of the denomina- 
tor, the critical coupling constant A* is fixed as a function 
of T, /i q , and B in unit of A, that is, 



1 A^c v-, „, /"da; v- / 2\q f B\n 

j n — U 



X 1 



i 



l 



e /3(a;A-^ q )_|_^ c f3(xA+n q ) ^± 

/3xAe /3(xA-^ 1 ) p xAe P(xA+^) 



(• e /3(xA- Mq ) + 1 )2 ^ e /3(xA+^) + iy 



(54) 



The numerical results are shown in Fig. [3] at zero den- 
sity (upper) and at finite density (lower) . We here like to 
point out that the numerical results are entirely smooth. 
This nice feature is attributed to the choice of the op- 
timized cutoff function. In fact, the Landau level sum- 
mation stops at n = Nk,B but the contribution in the 
vicinity of this upper bound is vanishingly small due to 
the weight y/l — 2\qfB\n/k' 2 . In other words, the opti- 
mized form (|32[) is a small regularization function, though 
it involves Heaviside's step function. 

In the upper panel of Fig. [3] we can confirm the Mag- 
netic Catalysis from the behavior of decreasing A* toward 
zero with increasing \eB\ at T — 0. In view of the upper 
panel of Fig. [3] however, the Magnetic Catalysis, A* = 0, 
is slowly reached for asymptotically large B only. We al- 
ready mentioned that the LLLA is exact for 2\qB\ > A 2 . 
In fact, the curve at T ~ in Fig. [3] is not the result 
at strict zero-T but at T = 10~ 5 A for numerical stabil- 
ity. As soon as finite T is introduced, A* is significantly 
pushed up and the Magnetic Catalysis can occur only in 
the limit of B —> oo, though A* = for any B at strict 
zero-T. 

It is also evident from the lower panel of Fig. [3] that 
the Magnetic Catalysis is lost immediately at finite /x q 
and A* never approaches zero even at T = then. More 
interestingly, we can see in the figure that A* increases, 
meaning that chiral symmetry breaking weakens, with 
increasing B contrary to the trend at fi q = 0. These 
are quite suggestive results implying that the modifica- 
tion in the QCD phase boundaries induced by B weakens 
at higher temperature, and B would rather favor chiral 
restoration for cold and dense quark matter. 
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FIG. 3. Critical coupling constant A, of chiral symmetry 
breaking and restoration as a function of the temperature T 
and the magnetic field \eB\ in unit of A for zero quark chem- 
ical potential fi q = (upper) and /i q = 0.3A (lower; with 
zero-density results shown for reference). 



IV. NON-LOCAL VERTEX AND 
RESUMMATION 



So far, we have worked in the leading order of the 
l/-/V c -expansion which amounts to neglecting the dia- 
grams (b),(c),(d) in Fig.[T] Once they are included in the 
calculation, the overall coefficient is modified from 2N C 
to 2N C + |, leaving aside the complexity due to isospin 
symmetry breaking by B. 

Hence, in practice, the sub-leading contributions have 
the minor effect of changing the overall factor. Strictly 
speaking, however, this is not the end of the story. The 
running coupling constant has an interpretation in terms 
of the meson propagation in the intermediate state as 
illustrated in the schematic representation (I3TJ1) . Because 
we defined A^ from the four-point vertex function with 
zero momentum insertion from external legs, we could 
describe the meson propagation only for zero momentum, 
as explicitly indicated in Eq. (|M|) . Such a treatment is 



correct for (a), but not adequate for (b), (c), and (d). 
In this section we will discuss qualitative effects of the 
momentum-dependent (i.e. non-local) vertices mediated 
by the meson propagation. 



Partial resummation 





Using the graphical representation (1301 we can express 
the diagram (d) as a box-type diagram as follows; 



(55) 



Here, it is clear that two meson propagators must carry 
the internal loop momentum and \k in the diagram 
should be replaced by the momentum-dependent cou- 
pling Afe(p). For the diagrams (b) and (c) as well as 
(a) we can utilize a similar representation with non-local 
vertices. 

One might have thought of a different type of the box 
diagram, that is, the one with the s- and i-channcls 
flipped. However, such a diagram does not give an ad- 
ditional contribution because, as we have explained be- 
low Eq. (j2Tj) . the color structure is decomposed into two 
pieces. We are now handling the part of 5 al3 5~ fS and the 
flipped diagram belongs to another part of S^S^ 1 . This 
is a crucial difference from the box diagram in terms of 
flavor indices encountered, for example, in the QED anal- 
ysis [i3 |. 

To take account of these effects properly, therefore, we 
should use the momentum-dependent coupling constant 
and it is natural to anticipate the following replacement; 



Afc 



A A 



i-A A [n fe (o) 

A fc (p) 



n A (o)] 
a a 



(56) 



i-A A [n fe (p)-n A (p)] 



This is nothing but the coupling arising from a full 
momentum-dependent bubble resummation and can be 
achieved within 2PI-approximations to the flow (4CJ] . 

Within such a momentum-dependent approximation 
we should discuss anew the validity of the LLLA. We 
can imagine the complexity from a simple fact that the 
energy dispersion of charged quarks is Landau quantized, 
while neutral mesons consisting of those quarks do not 
feel B at all and thus no Landau quantization occurs 
if they are tightly bound. We can indeed understand 
this from Eq. (l56l) ; the point is that Hk(p) in Eq. (|56l) is 
the polarization in the presence of the IR cutoff fc, that 
is, the IR modes with p 2 + 2\qfB\n < k 2 are massive 
and do not run in the loop, which makes Iik(p) depen- 
dent on k. As long as \qfB\ is well smaller than A 2 , 
the magnetic field cannot resolve the internal structure 
of bound-state system and the LLLA is not sufficient. In 
this case the apparent suppression of quark propagator 
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by \q/B\ is canceled by the energy denominator in the 
polarization calculation. On the other hand, if \q/B\ is 
greater than A 2 , all non-zero Landau modes (n ^ 0) are 
outside of the IR cutoff region and only the zero-mode 
can depend on k. Therefore, only the zero-mode survives 
in Ilfe(p) — IIa(p). Roughly speaking, the transverse mo- 
tion of neutral mesons gets frozen if the magnetic field 
strength is larger than the inverse squared of the meson 
size and this could be regarded as -B-induced dissociation 
or deconfinement. Such a transitional change of meson 
structure as a function of B is a very interesting theoret- 
ical question beyond the scope of this present work. 

Let us just assume that \q/B\ is sufficiently large and 
the LLLA holds in the computation of Hk(p) — ILv(p)- 
Then, after some calculations, we find, 



n fc (p) - n A (p) 



2tt 2 ^ 

/ 



W 



B\ 



\"Pz 



3pl 



ln ^ + 1f + 8§ 



(57) 



where we expanded the results in terms of p^/k and p z /k 
and took A ^> p 4 , p z to neglect P4/A and p z /A. We 
can immediately confirm that we can recover the run- 
ning coupling (|46)l by plugging Eq. (fSTf into Eq. (|56|) 
at zero momenta. Because p z is the loop momentum in 
the diagrams (b)-(d) which is cutoff at k by the quark 
propagator, p z and p 4 are certainly smaller than k. Since 
the introduction of k explicitly breaks Lorentz symmetry, 
Eq. ([57)1 is a function of not only p\ + p\ but also \p z \ 
independently. We have introduced r k (p) in the quark 
propagator to suppress the IR modes with \p\ < k and, 
somehow, the meson propagator should have the IR cut- 
off generated dynamically. Then, the spatial momentum 
of meson is shifted up by the IR cutoff k. To take account 
of this cutoff effect qualitatively, it would suffice to keep 
the linear term only in Eq. (|57[) . The running coupling 
with momentum dependence is then expressed as 



A fc (p) 



A» 



1 + 



N c X k 
Air 2 



f 



\Pz 



(58) 



We can then perform the momentum integration. For 
qualitative discussions we can treat X k (p) with some typ- 
ical p in Eq. ([58"]) . 

In the vicinity of chiral phase transition, A^ in \k(p) 
becomes larger and larger as k decreases, which is a flow 
pattern implied from Eq. (|45j) . Naively, one might have 
thought that the effect of the diagrams (b)-(d) is just 
to replace the overall factor in the /^-function, but this 
is not the case once the momentum dependence in the 
running coupling constant is taken into account. As is 
clear from Eq. ([5"5]) . we can see, 



Afc(p) -► (const.) (A fe 



0) 



(59) 



This means that the contribution from the diagram (a) 
which is proportional to A 2 overwhelms (b)-(d) which are 



proportional to Afc(j?)Afc and X k (p) 2 , respectively, near 
the chiral phase transition. Therefore, our analysis of (a) 
in the previous section is not changed by the inclusion of 
(b)-(d) qualitatively. 

To go into more quantitative studies, however, it is 
inadequate to use X k that is solved previously, but we 
need to solve the flow equation self-consistently including 
the non-local vertices. 



B. Self-consistent resummation 

Finally we advance our truncation scheme for taking 
account the non-local vertices X k (p)- The result for Afc(p) 
discussed previously stems from a standard resummation 
of bubble diagram in the s-channel or the leading-order 
contribution in the 1/N C expansion. Now we utilize the 
flow for going beyond this approximation. For that pur- 
pose we parameterize the full s-channel coupling X k (p) 
as 



Afc(p) 



A fc (p) 



i-A fc (p)An fc (p) 



(60) 



with an unknown function \k(p)- One can understand 
this Ansatz as a running coupling (1561) with shifted renor- 
malization. That is, one can immediately reach Eq. (|60[) 
by the following replacement in Eq. (I5t 



(61) 



-L + An A (p) -> * 

Afc \k{p) 



We note that we have to promote the right-hand side as 
a momentum-dependent quantity to take account of the 
full resummation including (b)-(d) in addition to (a). We 
can fix the normalization at vanishing momentum as 



A fe (0) = A fe (0) = A fe 



(62) 



The flow of Afc(p) is, in principle, given by the flow equa- 
tion with diagrams having momentum insertion from the 
external legs. For the derivation of the flow of X k (p) we 
can take the derivative on Eq. (jSU)) to find, 



9fcA fc (p) _ d k \ k (p) 
A&0) 



1 



A fe (;p) l + \ k ( P )ATl k (p) 



- X k (p)d k AU k (p) . 

(63) 

For the solution of the flow d k X k (p) in Eq. (|6"31 we have 
to insert the flow of dkX k (p) in a given approximation. A 
fully self-consistent flow results in loops with momentum- 
dependent four-point vertex which would require numer- 
ical computations of the Dyson-Schwinger type. 

In the present work we invoke further approximations 
motivated by the physics under discussion. This also 
facilitates the task of solving the flows and to uncover 
the underlying mechanisms. 

The coupling parameter Afc (p) encodes the pion effects 
beyond the leading s-channel resummation. These ef- 
fects, as is already the case in the previous subsection 
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for the s-channel resummation, are suppressed for non- 
zero internal momenta. This structure enables us to solve 
the flow for A& (p) within a derivative expansion in the s- 
channel momentum p. It also helps that the successive 
integration in terms of k already includes the momentum- 
dependence of Afe (p) within the identification, 



Xk=a(q) — Afc =(J (0) 



(64) 



Equation (|6"4"|) also provides a self-consistency check of 
the current derivative expansion. In the present work 
we take the lowest-order derivative expansion as A& (p) ~ 
Afe(O) = Afe [see Eq. ([62]) ] and, analogously to the flow 
equation (|23[) . we arrive at 



9feAfe 



E 

i—a,b.c.d 



Diag 



« 



(65) 



where Diag^ with i — a,b,c,d stands for the momen- 
tum integrations corresponding to the diagrams (a)-(d). 
With the coupling Xk(p) put down in Eq. ([60|) the dia- 
grams take the form; 



d 4 p 






X\ h(j?4,p) , 



(66) 



f=u,d' 



Diag (6) + Diag 

f=u,d 



(c) 



A 2 fe 



d 4 p 
(27) 1 l-A fc An fe (p) 



h(P4,P) , (67) 



-W-£/<& 



A? 



f=u,d' 



(2tt)4 [l-A fc An fc (p)]2 



h{pA-,p) 



with the common integrand given by 



h(P4,P) 



dkrkipjP^il + rkif)} 



{Pi 



■p[i 



r k (p)] 2 Y 



(68) 



(69) 



For the diagram (a) the s-channel momentum is the ex- 
ternal one for both couplings and hence A^ appears in 



the integration just in the same way as in Eq. (|23|) . For 
the diagrams (b) and (c) one of two couplings has zero 
momentum in the s-channel and the other has the loop 
momentum p. Finally, for the diagram (d), we have p for 
both vertices. 

It should be noted that we have neglected isospin sym- 
metry breaking in the above expressions. In other words, 
we here introduced simplification by setting the electric 
charges of w-quarks and <i-quarks to be equal; q u = qa- 
Without this, the polarization Hk(p) is dependent on 
which flavors of quarks each diagram involves and thus 
Xk(p) should be also flavor dependent. As a result, for 
example, the charged pion propagations are suppressed 
and the contributions from (b) and (c) are vanishing due 
to cancellation between the ir° and a processes. We could 
have discussed all these complexities here, but they are 
not the main subject of the present work. We are only in- 
terested in formulating the resummation procedure here 



-200 
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No Resummation \ 
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FIG. 4. Flow patterns of the four-Fermi coupling constant 
without resummation shown by the dotted line (in which the 



overall coefficient is 2N C 
by the solid line. 



+ 4) and with resummation shown 



and looking over qualitative effects on the flow. Though 
we have to leave from real electromagnetism for the mo- 
ment, we shall stick to simplicity with q u — qd- 

Now we have a closed flow equation for A^. It takes 
into account the scale dependence of the coupling at van- 
ishing momentum and goes beyond the standard bubble 
resummation. The current approximation fully restores 
RG-invariance; the coupling Xk{p) is invariant under a 
RG rescaling and looses all reference to the coupling at 
the cutoff, i.e., d\ A Xk=o = 0. 

Let us now discuss how the /3-function is modified by 
this prescription of resummation. For this purpose we 
will again focus on the calculation in the strong B limit 
which enables us to use the LLLA for analytical calcula- 
tions. In the LLLA, as we have seen, ALIfc(p) can be ap- 
proximated as AU k (p) ~ -{N c /^ 2 )Y, f \qfB\{\p z \/k). 
Since the momentum dependence is such simple, one can 
easily integrate each Diag^ with non-local vertices. 

After the ^-integration, with this resummation proce- 
dure, the /3-function changes from Eq. ([4"T)l into 



d k X 



(B) 



■kA t 



-^c(Af) 

X (B) 



2 A fc 



(B) 



(AD 5 



(70) 



81 



3 X (B) 
2 A fc 



The behavior of the right-hand side is shown in Fig. |U 
We can see that the resummation has a minor effect as 
stated in qualitative discussions in the previous subsec- 
tion. The resummation tends to disfavor chiral symmetry 
breaking slightly, which is consistent with the intuition 
that meson fluctuations would rather restore chiral sym- 
metry. 
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V. CONCLUSION 

We have investigated the RG flow pattern of the four- 
Fermi coupling constant A& and clarified the intuitive pic- 
ture of the chiral phase transition at zero and finite tem- 
perature and baryon density. We have applied the for- 
malism to study the effects of the homogeneous magnetic 
field. In the RG analysis it is extraordinarily simple to 
understand the phenomenon called the Magnetic Cataly- 
sis. If the magnetic field is strong enough, the transverse 
dynamics of charged particles is completely frozen and 
the dimensional reduction effectively occurs. Then, with 
the IR cutoff scale k, the transverse phase-space factor k 2 
is replaced by the Landau degeneracy \qfB\/(2%), which 
changes the flow pattern significantly. 

For dimensional reasons, in the absence of an external 
magnetic field B, a dimcnsionless coupling constant is 
given by Afcfc 2 , which tends to go to zero as k decreases 
unless the initial Aa exceeds a critical value. In contrast, 
if strong magnetic fields B are applied, the dimension- 
less combination of coupling constant is rather character- 
ized by (%2* \qfB\/2ir )\k- For this combination one can 
show that its RG flow leads to a diverging coupling how- 
ever small the initial Aa is. This is the simple RG-picture 
of the Magnetic Catalysis. 

Our quantitative studies in numerical calculations have 
made it clear that the dimensional reduction, however, is 
a sensible approximation only when B is unrealistically 
large and the temperature T and the chemical potential 
jj, q are sufficiently close to zero. The logarithmic singu- 
larity that is responsible for the Magnetic Catalysis is 
exactly canceled by another one from the matter parts 
at finite T and/or /x q . 

Finally, we have included the effects of non-local ver- 
tices mediated by the meson propagation with finite mo- 
mentum. The resummation procedure argued in this 
work thus deals with separate contributions from me- 



son loop effects and contains non-trivial contributions 
through the RG flow; see Ref . [4(j for a similar kind of re- 
summation. The related set of diagram has no qualitative 
impact on the location of the chiral phase transition. The 
present novel approximation to the NJL-model takes into 
account the back-reaction of intermediate mesonic states 
and hence relates to the full quark-meson flow in (P)QM 
models; see e.g. [121. [45144 7j . We hence expect the present 
approximation to be sensitive to the critical properties of 
the phase transition. The discussion of critical properties 
is deferred to future extensions of the present work. 

Although it was not quite necessary in the present 
work, it would be an intriguing question how to describe 
the neutral meson as charged quark composite in exter- 
nal B. Our RG analysis clearly shows that the lowest 
Landau level approximation works for the neutral me- 
son propagation only when B is greater than the cutoff 
scale A that roughly corresponds to the transverse me- 
son size inverse. This means that the quark dissocia- 
tion or deconfinement could be induced by the B effects, 
which would offer us a new opportunity to shed light on 
physics of quark confinement. Besides, it is an important 
but poorly understood problem to clarify the physical 
meaning of the conventional criterion for quark confine- 
ment using the Polyakov loop at strong _B, which would 
help us with understanding confinement/deconfinement 
at high density and hopefully more detailed characteris- 
tics of Quarkyonic Matter [48| . These are all interesting 
future extensions of the present work. 
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